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Jupiter’s upper atmosphere is composed of a neutral thermosphere and 
charged ionosphere. In the latter, the dominant molecular ion H,” emits 
inthe near-infrared, allowing for the remote exploration of the physical 
properties of the upper atmosphere. However, the Jovian low-latitude 
ionosphere remains largely unexplored because H,* emissions from this 
region are faint and spectrally entangled with bright neutral species, such 
as CH,. Here, we present James Webb Space Telescope H,* observations 

of Jupiter’s low-latitude ionosphere in the region of the Great Red Spot, 
showing unexpected small-scale intensity features suchas arcs, bands and 
spots. Our observations may imply that the low-latitude ionosphere of 
Jupiter is strongly coupled to the lower atmosphere via gravity waves that 
superimpose to produce this complex and intricate morphology. 


The extent to which the non-auroral upper atmosphere of Jupiter is cou- 
pled to the convective-radiative atmosphere below and to the magne- 
tosphere above is largely unknown because ground-based observations 
donot offer the sensitivity required to investigate this. The James Webb 
Space Telescope (JWST) performed observations of Jupiter’s Great Red 
Spot (GRS) with the NIRSpec Integral Field Unit (IFU) on 27 July 2022 as 
part of the Early Release Science programme (ERS #1373 (ref. 1)). The 
F290LP/G395H filter/grating produced a spectrum between 2.87 and 
5.27 um with a spectral resolving power of R = 2,700. In this Article, we 
use the short-wavelength (NRS1) detector only, covering 2.85-4.01 um. 
The observational setup is shown in Fig. la, with JWST acquiring six 
IFU spectral cubes over a period of 53 min, each with an effective inte- 
gration time of 5.4 min. The observational details are summarized in 
Table 1. The NIRSpec IFU has 30 x 30 spaxels, each measuring 0.1’ x 0.1’, 


producing a spatial resolution of ~300 km per spaxel on the planet. This 
wavelength region contains emissions from the molecular ion H,;",a 
dominant ion in the charged-particle ionosphere, along with H*, situ- 
ated at altitudes >300 km above the cloud deck’. This ion is produced 
by both electron-impact ionization about the magnetic poles and solar 
photo-ionization on the dayside of the planet, which dominates as a 
production mechanism in the low-latitude ionosphere’. In the absence 
of dynamic interactions with the atmosphere and magnetosphere, the 
low-latitude ionosphere may be expected to emit in a relatively uni- 
form manner, given that the incoming solar flux at these heliospheric 
distances will be constant across the disk at any given moment, witha 
peak production of H,” occurring at the subsolar point. At low latitudes, 
the H;* emissions are roughly ten times weaker than the auroral region* 
and are blended with strong solar excited non-local thermodynamic 
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Fig. 1|JWST observations of Jupiter’s GRS. a, The JWST NIRSpec IFU mosaic 
sequence used to capture Jupiter’s GRS area on 27 July 2022, shown together 
with the 3 pm aerosol scattering from the upper troposphere. The times are the 
midpoint of each observation, and the background image is a JWST Near Infrared 
Camera (NIRCAM) observation from 27 July 2022 20:27:10 using the F440W 
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filter. Figure 2 shows spectra at points P, and P,. b, The integrated H;* radiance 
extracted from the NIRSpec observations. The horizontal grey lines show the 
magnetic dip angle, and the vertical grey lines shows the orientation of the 
magnetic field lines as mapped above the surface, spaced by 10° longitude’. The 
oval shows the location of the GRS. 


equilibrium (LTE) methane emissions’. This makes any observations 
of this region from the ground challenging‘ and therefore a suitable 
target for JWST. The integrated ionospheric H,* emissions are shown 
in Fig. 1b, and the analysis procedure is outlined below. 

In Fig. la, two points are indicated: P, at the GRS and P, away from 
it. The NIRSpec spectra for these points are shown in Fig. 2b, together 
witha model H,* spectrum in Fig. 2a. The observed radiances are clearly 
much higher at the GRS than away from it, showing an enhancement 
in the reflective aerosols in the upper troposphere and lower strato- 
sphere at the GRS, extending up to 40-50 mbar (ref. 7). This figure 
also illustrates the challenge in extracting the H,* spectrum from this 
region, which is rich in spectral features: the H,” signal is very weak 
compared to bright structures generated in the lower part of the atmos- 
phere, with H;* peaking at -1 mW m”? sr‘ wm ‘and the peak radiance at 
3.9 um reaching up to -20 mW m” sr? pm“ where theH;" Q(1, 07) line is 
situated. However, between 3.2 and 3.6 um, deep CH, absorbs most of 
the incoming sunlight, rendering this region spectrally dark, and the 
discrete H,° emission can be readily extracted after the subtraction of 
the high-altitude solar-pumped non-LTE CH, spectrum. This process is 
repeated for each spaxel and is illustrated in Fig. 3 and expanded upon 
in the Methods. Note that the strong H," line at -3.46 um is omitted in 
this analysis because it is blended with a bright feature produced by 
the absence of methane opacity in the lower atmosphere’, whichis not 
readily modelled. The spectral regions used in this study are shaded 
grey in Fig. 2, containing 21H,” emission lines (Supplementary Table 1). 
Overall, the solar-pumped CH, emissions do not show any substantial 
small-scale variation across the field of view (Supplementary Fig. 1). 

The H,* emissions extracted from the NIRSpec data can be seen in 
Fig. 1b, together with the magnetic field topology across this region’. 
Astriking result of JWST observations that was never seen in previous 
ground-based images’ is the presence of intricate small-scale struc- 
tures. These are visible in this small patch on the planet, covering only 
37° longitude and 27° latitude. In fact, ground-based observations 
need to average over several nights to bring out low-latitude features, 
preventing us from seeing the instantaneous ionospheric structure, 
which now reveals its complex morphology. There are discrete fea- 
tures that form bands and arcs ~1,000 km wide and up to ~6,000 km 
long. On average, H,’ is relatively dark inside the GRS, with two bright 
regions to the northwest and southeast of it. Across the field of view, 
the integrated radiance varies by up to a factor 4. The intensity of 
H;* is governed both by ion density and by the temperature of the 
upper atmosphere, with which the ion thermalizes. Via H," spectral 
retrievals, these two parameters can be disentangled and the driver of 
the ionospheric emission features can be determined. The extracted 


Table 1| Summary of the JWST NIRSpec IFU observations 


Nbr Date start Time end Time mid Lon Lat LT (h) 
(UTC) 

1 27 July 2022 09:12:48 09:10:02 303.8 -177 11.0 
09:07:15 

2 27 July 2022 09:22:17 09:19:31 298.1 -11.6 11.8 
09:16:44 

3 27 July 2022 09:31:46 09:29:00 290.8 -14.6 12.7 
09:26:13 

4 27 July 2022 09:41:15 09:38:29 295.6 -20.7 12.7 
09:35:42 

5 27 July 2022 09:50:44 09:47:58 2871 -241 13.7 
09:45:11 

6 27 July 2022 10:00:13 09:57:27 282.5 -17.9 14.4 


09:54:40 


Tile acquisition times at Jupiter's GRS and the tile centre planetocentric latitude (Lat) and 
System III west longitude (Lon), as well as the local time (LT). The observing sequence is 
shown in Fig. 1. 


spectra for each spaxel (for example, Fig. 3c) are fitted using the open 
source H, Python package h3ppy, and the results of these retrievals 
are shownin Fig. 4. The median temperature of all spaxels is 713 + 20 K, 
with a median density of 4.5 + 0.4 x 10 m”, broadly consistent with 
previous studies™®™”. 

In Fig. 4, the oval indicates the location of the GRS (from Fig. 1), and 
the contours show a Gaussian smoothing with a width of 2° longitude 
and latitude of each fitted parameter, showing the broad-scale struc- 
ture present inthe data. The temperature distribution in Fig. 4a shows 
higher temperature toward higher latitudes, and three regions appear 
hotter, with a median temperature of 738 K: (1) at the centre of the GRS 
and (2) northwest and (3) southeast of it. Northwest of the GRS, the 
troposphere is extremely turbulent” and could be a source of gravity 
waves. The elevated temperatures at the GRS have larger uncertain- 
ties associated with them due to the elevated background emissions 
(Fig. 4c), but this could be related to an increased homopause altitude 
above the GRS, which destroys colder H,* at low altitudes via charge 
exchange, producing column observations of high temperatures and 
low densities, matching what is seen in Fig. 4. The higher noise levels 
observed in the GRS might be associated with the presence of elevated 
hazes that increase the brightness of the GRS at 3.3-3.6 um. The H," 
column-density retrievals in Fig. 4b show patchy structures that are 
largely mirrored in the integrated radiance (Fig. 1b). Next, we quan- 
titatively investigate the drivers of the discrete H,* radiance features. 
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Fig. 2 | JWST NIRSpec Jupiter spectra. a, A model H, spectrum generated with 
values consistent with low-latitude emissions (T=700 K and N=5 x10" m?). b, The 
NIRSpec IFU F290LP/G395H NRS1 spectra from individual spaxel points P and P, 


shown in Fig. 1, from inside and outside the GRS. The vertical green lines indicate 
the locations of the brightest H,* lines. The grey shaded area shows the wavelength 
region considered for this analysis, from which the H;” spectrum is extracted. 
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Fig. 3 | H,* extraction procedure. The analysis sequence applied to retrieve 
H;* model fits to the NIRSpec IFU data. a, The spectrum from point P, in Fig. 1, 
from which regions with bright H,* emission lines are extracted. b, ACH, model 
is fitted to the extracted spectrum. c, Once the CH, model is subtracted from 
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the data, isolating the H;* emission lines, the spectrum can be fitted using the 
Python package h3ppy to retrieve the temperature of the upper atmosphere and 
ionospheric H;* density. 


From the set of NIRSpec observations presented here, mosaic 
tile (vi) (Fig. 1) contains the clearest (dark) banded structures and is 
located away from the region with large errors at the centre of the 
GRS and regions of localized heating. To investigate whether these 
are produced by reduced temperatures or lower H,* densities, the 
relationship between these parameters and the integrated H,* radi- 
ance are plotted in Fig. 5. A third-order polynomial is fitted to each 
correlation (black line), with the shaded region showing the extent 
of the median uncertainty. The relationship between temperature 
and radiance in Fig. 5b shows that the temperatures are not strongly 
correlated with the observed H,° radiance, having a Pearson correla- 
tion coefficient of 0.15. However, Fig. 5c shows that as the radiance 
increases, the column density also increases, well outside the density 
uncertainties, with a correlation coefficient of 0.66. This shows that 
the observed changes in the H,* radiance is driven primarily by the 
H,* column density. The correlation relations for the other six mosaic 
tiles are shown in Supplementary Fig. 2. The Pearson correlation coef- 
ficients indicate that the radiance is better correlated with density 
than the temperature for all but one tile (tile (iii)). This may suggest 
that although the radiance is generally driven by density, there are 
small areas where the temperature is a significant driver, perhaps 
driven by gravity-wave heating”. 

Figure 5 also shows that higher temperatures are generally asso- 
ciated with lower column densities. This anticorrelation is well estab- 
lished and expected (for example, refs. 14,15) and is produced by the 
fact that the intensity of H,° emissions is driven exponentially by tem- 
perature and linearly by density”, so they are intrinsically anticorre- 
lated. In general, the spatial differences in temperature and density can 
be trusted if the differences are larger than the derived uncertainties“, 


which is not the case for temperatures in Fig. 5b but is the case for 
densities in Fig. 5c. 

In Fig. 4b, the H,* column densities vary by a factor of ~4. In the 
polar region, the column densities are higher by a factor of ~10, but 
the column-density variability across the H,* aurorae varies by similar 
factors”. The H,* density variations at high latitudes are related to the 
discrete auroral field-aligned currents generating vast amounts of 
H," via impact ionization. Such large density gradients are a surprise 
at low latitudes, however, where ion production is dominated by the 
smoothly varying incident solar flux. High-altitude in situJuno plasma 
observations have also suggested that there are significant fluctuations 
in the H;* density at low latitudes’®, which may be related to the density 
variations seen here. 

The results of this study form the starting point for understand- 
ing the short-term and small-scale drivers of Jupiter’s low-latitude 
ionosphere. At Earth, the low-latitude ionosphere is dynamic and rich 
in features, dominated by electrodynamics”. Neutral winds play an 
important role inthe behaviour of this region, carrying the ionospheric 
plasma with them away fromthe subsolar point. This generates a strong 
duskward electric field (E), known as the equatorial electrojet. Close 
to the equator, the magnetic field (B) runs parallel to the surface and 
the magnetic field-line dip angle is zero, also called the magnetic dip 
equator. At this location, the E x B force moves ions up in altitude, gen- 
erating the ‘plasma fountain’, producing a depletion of plasma at the 
magnetic dip equator and enhancements in plasma density at either 
side of it as the plasma diffuses back down due to gravity, known as 
the equatorial ionization anomaly (EIA). The magnetic field of Jupiter 
has a magnetic field polarity opposite that of the Earth, reversing the 
direction of B at the equator. This also likely reverses the equatorial 
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Fig. 4 | Temperature and H;* ion density retrievals. The physical properties 
ofthe ionosphere >300 km (ref. 2) above the Great Red Spot area. a, the 
temperature, T, of the upper atmosphere, b, H;* column integrated density, N, 
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c, the uncertainty on the retrieved temperature, AT, d, the uncertainty onthe 
retrieved H;* column integrated densities, AN. a and b also show contours of the 
smoothed parameters to reveal broader-scale structures. 


electrojet to run dawnward, still resulting in an Earth-like upward E x B 
force, producing a plasma fountain, removing H,’ ions from the mag- 
netic dip equator. However, our observations show no clear decrease 
in the H," density close to the dip equator, nor any evidence for EIAs 
(Fig. 4b and Supplementary Fig. 3). 

Another important driver of ionospheric density perturbations 
at Earth are gravity waves, which produce travelling ionospheric dis- 
turbances (TIDs), where a single travelling wave produces ordered 
density perturbations in the ionosphere”. Above Jupiter’s GRS, the 
ionosphere is likely subject to a complex barrage of upward travelling 
gravity waves originating from different tropospheric storms” ~*. 
This process might explain the complex structures see in Fig. Ib—the 
superposition of multiple gravity waves changing the H," density of the 
ionosphere—and would also explain the arc-like structures seen in the 
data. For example, the right-angled feature in mosaic tile (vi) could be 
the result of two perpendicular waves. 

A further process that shapes the ionosphere of the Earth at lower 
latitudes is the propagation of Rayleigh-Taylor instabilities, espe- 
cially in the post-dusk sector”. This process is capable of generating 
extended features and requires a trigger, such as the removal of the 
low-altitude ionosphere after dusk or a density fluctuation caused 
by gravity waves. However, the main H;* density peak at Jupiter is in 
photochemical equilibrium”, which will act to limit the development 
of Rayleigh-Taylor instabilities on the dayside; the relatively short 
lifetime of H,° means that these cannot be propagated effectively. 
Therefore, we suggest that TIDs generated by gravity waves are the 
most likely mechanism responsible for the small-scale features seen 
inthe JWST data. 

At Earth, another source of ionospheric irregularities is meteoritic 
seeding of the upper atmosphere that produces vertical stratification 
of the ionosphere, knownas sporadic E layers, which generates local- 
ized increases in the ionospheric density”. Jupiter’s large mass makes 
these impacts more common than at the Earth and could therefore 
be an important contributor to shaping the ionospheric structure, 
manifesting as H;* density perturbations close to the homopause”. 
The meteoritic ions (Mg* and Fe’) are long lived, with lifetimes of ~300 
Jovian rotations; the sporadic layers can also be shaped by gravity 


waves, and the spatial variability of electron density can significantly 
alter the lifetimes of H,”, shaping the observed densities. 

Gravity waves in Jupiter’s atmosphere, analogous to TIDs at the 
Earth, have been modelled”, showing that they can produce H;* den- 
sity variations of over an order of magnitude, as well as much smaller 
perturbations in temperature, up to ~30 K. Because the temperature 
uncertainties in Fig. 4c are of a similar magnitude, we expect to resolve 
gravity-wave perturbations only in the H;° column density and not in 
temperature. The dominant ions in the ionosphere of Jupiter are H* 
and H,", with lifetimes of ~1Jovian day and ~1,000 s, respectively, at 
~300 km (ref. 21). For waves with periods close to the H; lifetime, only 
very small H;* density perturbations are expected; but for longer waves 
with periods closer to the H+ lifetime, significant H,* variations are pre- 
dicted. This is because the chemical reactions that govern H,_ (thatis, 
production and destruction) occur on much shorter time scales, which 
means that H,° remains in chemical equilibrium along the wave, so the 
H,’ radiance features are a result of the gravity-wave-induced density 
perturbations in the H, source population. The gravity waves seen in 
the Galileo entry probe thermal profile” are predicted to generate 
H,’ density variations of 5% to 23%, whereas here we see variations of 
up to a factor of ~4 (for example, see Fig. 5c). Note, however, that the 
modelled H;* density variations are highly dependent on the wave 
parameters used, and they do not consider the superposition of mul- 
tiple gravity waves. 

The temporally averaged view of Jupiter’s ionosphere produced 
by adding almost 50 nights of observations* showed that the intensity 
of H;* is strongly linked to the topology of the magnetic field, suggest- 
ing that electrodynamic processes shape the structure of Jupiter’s 
low-latitude ionosphere. In stark contrast, the JWST observations 
presented here suggest that the instantaneous structure of the iono- 
sphere is largely driven by wave dynamics. The equator-ward edge 
of IFU mosaic tile (ii) reaches the magnetic dip equator (Fig. 1b), but 
there is no H," density decrease present to suggest the existence of a 
plasma fountain or an EIA (Fig. 4 and Supplementary Fig. 3). This is 
not, however, inconsistent with the electrodynamics revealed by the 
long-term study. Wave dynamics must govern the short-term behaviour 
of the low-latitude ionosphere, and only when these are averaged out 
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Fig. 5 | Drivers of H;* radiances. a, The integrated H,* radiance for NIRSpec IFU mosaic tile (vi) (Fig. 1). b,c, Correlation between temperature and radiance (b) and 
column-integrated H,* density and radiance (c). The black line shows the third-order polynomial fit to the data, and p indicates the Pearson correlation coefficient. 


over longer time scales can we reveal features associated with electro- 
dynamics. In other words, in the low-latitude ionosphere of Jupiter, 
wave dynamics appear to dominate over electrodynamics inthe short 
term, at least in the GRS region. The opposite is true at Earth”. Addition- 
ally, Jupiter’s low-latitude ionosphere has long been considered quite 
bland”, especially when contrasted to the dynamic auroral zones. The 
observations presented here show that this is far from true, and there 
is a richness in features that have never been seen before. 

Although solar photo-ionization is expected to dominate the 
production of H;* at low latitudes’, there has been speculation that 
particle precipitation could also occur in this region”. An increase in 
particle impact ionization would serve to increase the observed H," 
density, due to increased ion production. However, the very complex 
small-scale pattern seen in Fig. 1b would imply an equally complex set 
of field-aligned currents across this small region, and there is hitherto 
no evidence of such currents away from the auroral region. 

One of the largest outstanding questions concerning the upper 
atmospheres of the giant planets is the energy crisis, wherein the 
observed low-latitude temperatures of these regions are hundreds of 
kelvin hotter than what the solar energy flux can provide (for exam- 
ple, ref. 28). The two main candidate resolutions to this crisis are 
equator-ward redistribution of auroral heating and the breaking of 
acoustic or gravity waves generated in the turbulent troposphere. Both 
of these potential solutions have challenges to overcome. According 
to some models”, the movement of energy from the pole down to the 
equator is largely prohibited by the strong Coriolis forces generated by 
the large size and rapid rotation of these planets, and the effectiveness 
of gravity waves as a heat source remains largely unknown. Pole-to-pole 
observations of the upper atmosphere of Jupiter and Saturn”? show 
that both have a temperature gradient that very gradually decreases 
fromthe pole to the equator, indicative of global transport of heat away 
from the auroral region. The somewhat higher temperatures seen at 
higher latitudes in this study are broadly consistent with this scenario. 

At Saturn, momentum loading of the upper atmosphere, which 
could be driven by gravity waves, has been proposed as a mechanism 
that could reduce the efficacy of the Coriolis force in preventing the 
movement of heat from the pole down toward the equator”. The dis- 
turbances observed here are highly suggestive of gravity waves present 
in Jupiter’s ionosphere. Having an upper atmosphere replete with 
these waves could similarly reduce the efficacy of Coriolis forces at 
Jupiter, enabling the equator-ward transport of heat, explaining the 
gradual pole-to-equator temperature gradient reported previously", 
suggestive of a global redistribution of auroral heating. In addition, 
the large-scale wave of elevated temperatures previously seen at the 
polar region” may also contribute to the wave activity seen at lower 
latitudes when it propagates towards the equator, although this was 
largely affecting the observed temperatures, not the densities, as is 


seen here. Although there is observational evidence of redistribution 
of auroral energy, reproducing these results with numerical model- 
ling is very challenging, and the exact mechanism that allows heat 
to be redistributed is still an open question. By varying the incoming 
solar-wind pressure, and therefore the auroral heating, modelling has 
been able to show that some energy is able to move equator-ward”’, 
albeit only at levels of tens of kelvins. 

The temperatures derived from the JWST data are relatively 
high, as expected from the ‘energy crisis’, but they are also very uni- 
form across the field of view. Although heating has previously been 
observed above the GRS”, no extreme temperatures were found inthe 
present study. This potentially indicates that extreme GRS heating is 
time-variable or that upwelling from below can occasionally carve out 
the bottomside of the ionosphere”, resulting in only high-altitude, 
high-temperature H,* remaining. However, the presence of multiple 
interfering gravity waves in the JWST data may have important implica- 
tions for ionospheric heating. Turbulence is present in the troposphere 
almost everywhere, so the breaking of gravity waves could produce 
localized heating at a large number of locations across the planet. As it 
stands, there is no modelling consensus to what degree gravity waves 
heat the upper atmosphere”. Further detailed modelling is required to 
determine the physical properties of the gravity waves responsible for 
shaping the low-latitude ionosphere: only certain combinations of wave 
parameters (for example, amplitude and wavelength) can explain the 
complex structures seen here, so these observations impose important 
constraints on future modelling efforts. Ultimately, by deriving those 
wave parameters based on observations such as these, better estimates 
of wave-induced heating can be made. 

Inthe auroral region, the large-scale structure of H;° emissions is 
governed by the interaction between the magnetosphere, ionosphere 
and atmosphere. However, it is possible that gravity waves would still 
provide small-scale ionospheric disturbances at high latitudes that 
are superimposed on the auroral emissions, so that the polar iono- 
sphere is shaped by forcing both from above (magnetosphere) and 
from below (gravity waves generated in the troposphere). Overall, 
however, the polar ionosphere is predominantly shaped by magne- 
tospheric processes, whereas the low-latitude ionosphere is driven 
largely by the coupling with the lower atmosphere. This highlights 
the importance of the ionosphere as an interface layer in a highly 
coupled neutral atmosphere-charged-particle ionosphere—mag- 
netosphere system. 

The fact that the lower and upper atmospheres are so strongly 
coupled at Jupiter may also imply strong coupling at other giant plan- 
ets, with momentum transfer via gravity waves being an important 
mechanisms to shape the global ionosphere. Because these waves alter 
the ionospheric density, they also change the conductivity, which in 
turn will regulate the efficacy of any interaction with magnetospheric 
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currents, particularly at high latitudes. These JWST observations form 
the proof of concept for future investigations of this region: for exam- 
ple, if these ionospheric irregularities were monitored over time, the 
progression of these waves could be directly measured, providing a 
view of temporal variability and an opportunity to characterize the 
wave amplitudes and wavelengths. 


Methods 

The GRS of Jupiter was observed by JWST NIRSpec* on 27 July 2022 
using the IFU and the F290LP/G329H grating/filter as part of ERS pro- 
gramme #1373: ‘ERS Observations of the Jovian System as a Demonstra- 
tion of JWST’s Capabilities for Solar System Science’ (Co-Principal 
Investigators: I. de Pater, T. Fouchet). The IFU produces a 
three-dimensional cube with two spatial dimensions covering 3’ x 3” 
field of view, with each spaxel covering 0.1 x O. 1’ and a third spectral 
dimension covering 2.87-4.01 um (using the NRS1 detector only) ata 
spectral resolution, at wavelength A, of R = 2 x 2,700. The NIRSpec 
IFU data was calibrated using the Python jwst calibration pipeline 
(v.1.9.5.dev3+g22acc46) up to Level 2 with Calibration References Data 
System context defined by jwst_1045.pmap. The MAST observation ID 
is jw01373-0005 t004_nirspec_g395h-f290Ip. Geometry, for example 
latitude and longitude, was assigned to the data using the Python class 
JWSTSolarSystemPointing (https://github.com/JWSTGiantPlanets/ 
NIRSPEC-Toolkit). 

Six IFU mosaic tiles were acquired to capture the full extent of the 
GRS, illustrated in Fig. 1 with details listed in Table 1. Each tile is a single 
observation made up of 30 groups giving an effective exposure time 
of 322 s (5.4 min). The use of a large number of groups was specifically 
to capture the relatively weak H,” emissions witha high signal-to-noise 
(S/N). Between 3.3 and 3.6 um, the mean S/N is 55, whereas at 3.0 umit 
is 274 (Fig. 2). The full2 x 3 mosaic took a total of 53.0 min to complete, 
with a mean time separation of 3.9 min between each tile. Table 1 also 
lists the local time at the centre of each tile, showing that the prenoon 
to mid-afternoon sector was sampled. The maps in Fig. 1 and Fig. 4 were 
projected with a latitude and longitude resolution of 0.05°. 

The Level 2 calibrated NIRSpec data contains a considerable 
amount of arc-like fringing, which appears to be especially apparent 
where the radiance is low, which is the case for the regions that contain 
emissions from H,’. To address this, a flat field was constructed fromthe 
six mosaic tiles, with the flat field of each wavelength bin being equal 
to the median of the six mosaic tiles. This works well when the spatial 
variability across the bin is small, for example for the H,” bins, but does 
not workat all when there are bright discrete features, for example dis- 
crete cloud features at 3 um. The flat-field procedure for a wavelength 
bin that contain strong H,* emissions is illustrated in Supplementary 
Fig. 4. Ultimately, the process used to remove the fringing does clearly 
improve the fidelity of the data and does not introduce any significant 
artificial features in the resulting maps. 

Figure 2a shows a model H;* spectrum generated using h3ppy 
with parameters broadly appropriate for Jupiter’s low-latitude 
ionosphere (temperature T= 700 K, column density N = 5 x 105 m°? 
(ref. 11). Figure 2b shows the full NRS1 spectra from both outside and 
inside the GRS (points P, and P, in Fig. 1). A significant number of H,” 
lines are available in this wavelength range, but disentangling weak 
lines froma very bright background is extremely difficult. However, 
between 3.2 and 3.6 pm there is significant absorption of sunlight by 
CH, in the troposphere and stratosphere, rendering this region very 
dark (indicated as grey in Fig. 2b and shown in Fig. 3a). Within this 
band, there are two dominant spectral features: (1) non-LTE solar 
excited CH, emission contained within 3.2 and 3.4 um and (2) discrete 
emission lines from H,”. To isolate the H,* spectrum, the solar excited 
CH, spectrum needs to be subtracted. The full radiative transfer mod- 
elling required to model the full CH, spectrum’ is outside the scope of 
this study, and a simplified, yet effective, approach is adopted here. 
The spectral shape of the non-LTE CH, emissions can be approximated 


by modelling two separate LTE components: the fundamental transi- 
tions (v, > ground) and the hot-band transitions (v; + v, > v4). Using 
the HITRAN® CH, line list, the radiance of the two components can 
be modelled separately at a fixed temperature (300 K), as well as a 
third-order polynomial background, which works well over limited 
wavelength ranges. Here, we use four separate narrow wavelength 
windows for the CH, fitting: (1) 3.23-3.31 um, (2) 3.37-3.40 um, (3) 
3.40-3.43 um and (4) 3.528-3.56 pm. An example CH, model fit for 
point P, is shown in Fig. 3b. The morphology of the integrated CH, 
radiance can be seen in Supplementary Fig. 1. The residual between 
the NIRSpec data and the CH, model then contains the H,* spectrum 
(for example, Fig. 3c), which is then fitted with h3ppy. The fitting 
program implements the H,' line list of ref. 36 and the partition func- 
tion of ref. 37. 

The analysis sequence illustrated in Fig. 3 was repeated for all IFU 
spaxels in this dataset, generating 6,927 individual H," spectral fits. 
Note that this is larger than the nominal 5,400 number of spaxels con- 
tained within six NIRSpec IFU observations and is due to the fact that 
the JWST pipeline drizzles the cubes onto a grid that is slightly larger 
than the 30 x 30 spaxel IFU when using the ifualign coordinate system. 

The magnetic field modelling and field-line tracing for Fig. 1 were 
performed using the JupiterMag model”, using the latest magnetic 
field model derived from the Juno magnetometer observations from 
33 perijoves of the primary mission’. 

For a broader context, Supplementary Fig. 5 shows the magnetic 
field strength and the smoothed temperature distribution derived from 
the data, and Supplementary Fig. 6 shows the wavelength-integrated 
H,* total emissions” on top of the H,* imaging observations of ref. 4. 
Supplementary Fig. 2 shows the unprojected integrated H,” radiance 
and the correlation between the temperature and H, density for tiles 
(i)-(v). 

Supplementary Fig. 7 cuts along -18° latitude and 291°W longitude 
of the H,* radiance, temperature and column density, emphasizing 
that the radiance features are largely driven by column density, not 
by temperature. The small-scale features are clearly seen, with large 
variability present within these profiles. 

The gaussian filter applied to Fig. 4 and Supplementary Figs. 5 and 
6is the ndimage.gaussian_filter function available in scipy. 

The JWST ERS Jupiter programme also included both NIRSpec 
FL170P/G235H filter/grating and MIRI Channel 1 Short observations 
of the GRS, both capturing a number of H,° emission lines at ~2 um 
and ~5 um, respectively. These could have provided a means through 
which to examine the longer-term spatial variability of the ionospheric 
features, because they were obtained at different epochs. However, 
the large magnitude and intricate structure of the underlying signal 
(deep thermal emissions and aerosol scattering) and the low S/N of the 
H,* emission lines renders this mapping impossible. Only the NIRSpec 
FL290LP/G395H filter/grating was optimized for the relatively weak 
H," emissions and used an effective exposure time roughly ten times 
longer than the other settings. 


Data availability 

JWST data used in this study were obtained from the Mikulski Archive 
for Space Telescopes at the Space Telescope Science Institute (https:// 
archive.stsci.edu/), which is operated by the Association of Universities 
for Research in Astronomy, Inc., under NASA contract NAS 5-03127 for 
the JWST.JWST NIRSpec ERS #1373 observations of Jupiter’s GRS used 
here are available via MAST at https://doi.org/10.17909/yedm-f607. 


Code availability 

The open source Python H,’ fitting code h3ppy is available at https:// 
github.com/henrikmelin/h3ppy. The JupiterMag field-line tracing 
code is available at https://github.com/mattkjames7/JupiterMag. All 
remaining code used for this analysis of the JWST data can be found at 
https://github.com/henrikmelin/jupiter-grs-h3p. 
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